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Zero Of a Solution of a Differential Equation 

(Using Library Routine Fl) (SADOI or DOl) 

Closed, with one program parameter (see below) 

16 

0, 1, 2, 5 

S3-S5 

Program parameter for routine Fl (SEE Fl) 

Steps proportional to N(g). 

Routine Fl begins at h. 

Entry 



p+1 



50 eF 
50 pF 



[e) is to become zero. 



26 xF This routine begins at x. 



This routine enters Code Fl repeatedly > adjusting the size 
of the steps, until one of the vsuriabies y vanishes. 
The variable which is to become zero is selected by the 
program parameter « 

The auxiliary routine used by Library Routine Fl must be 
so designed that the scaled derivatives 2\y. ' are all 
proportional to a scale factor stored in go This scale 
factor is altered during the convergence to the zero, and 
must be restored if further steps are to be carried out. 
|N(g) I is never increased by this routine, so that if the 
steps are initially sufficiently small for accurate 
integration they will remain so. This means that several 
steps may be necessary to reach the neighborhood of the 
zero; when this has occurred, about six steps are usually 
sufficient for convergence o N(e) is normally reduced to 
or to + 2"^^$ if hy' > N(g) the error may be as large 



- 2 - 

as 2"^^ hy7N(g) but this would be abnormal. 
Integration proceeds in the direction in which |N(e)| 
decreases. If a non-zero minimum of |N(e)| is encountered 
control may be incorrectly returned to the main program, 
or the machine may get into a loop. This routine should 
therefore only be entered when the integration is near 
the required zero. The scaled derivative of the variable 
which is to vanish must be non-zero when this routine is 
entered • 
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k2 5L 
lA 59B5 




Set link etc. 
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U6 5L 










26 5L 




Waste 
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L5 F 


By 2 








iwD F 


From Ik^ 


2%-' to 
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L3 F 
10 S3 




Form 2"^^ - hy' 
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lA 15L 










32 F 


By 1 


Links Test for convergence 
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50 F 
LO 15L 


By 0' 


Test if y > hy' 
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S6 F 
32 lOL 
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S5 F 










00 S3 




Near zeros Use Newton-Raphyson 
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66 F 
79 Si+ 
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26 12L 










71 F 


From 7' 


Test whether to reverse 
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56 13L 
U S4 




Reverse step 
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ho Sk 
26 13L 


From 10 


Waste 




13 


00 S5 
50 15L 


From 11 


Enter Program Fl 
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26 S5 
26 5L 
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